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AN  EFFECTIVENESS  PREDICTION  MODEL 
FOR  SUBMARINE  MISSILE  SYSTEMS 

I.  INTRODUCTION 

One  measure  of  the  effectiveness  of  a  submarine  missile  attack  is 
given  by  the  radial  distance  between  the  target  and  the  missile 
detonation  point.  Varying  degrees  of  damage  will  be  inflicted  by 
simply  placing  the  missile  in  proximity  to  the  target  at  time  of 
detonation.  Accordingly,  range  from  weapon  to  target  is  perhaps  the 
most  important  single  parameter  affecting  missile  performance  in  any 
tactical  situation.  Under  realistic  operating  conditions, 
discrepancies  between  the  actual  and  estimated  values  of  range  will 
invariably  occur  because  of  missile  aimpoint  delivery  errors,  as  well 
as  errors  inherent  in  the  target  localization  process.  To  reliably 
predict  weapon  performance,  these  errors  must  be  accurately  modelled. 

This  report  describes  a  real-time  computational  procedure  for 
predicting  the  effectiveness  bf  a  submarine-launched  missile  attack 
against  prescribed  target  threats.  A  functional  relation  is  developed 
which  expresses  the  probability  of  target  damage  in  terms  of  missile 
aimpoint  delivery  error  statistics,  weapon  damage  characteristics,  and 
target  localization  error  statistics.  This  functional  relation  is 
subsequently  transformed  into  a  digital  computer  algorithm  suitable 
for  making  on-line  weapon  performance  predictions. 

In  section  II  the  problem  is  formulated,  and  equations  for  the 
probability  of  target  damage  are  developed.  Section  III  contains 
approximations  to  the  probability  distributions,  while  section  IV 
presents  closed-form  expression  for  probability  of  target  damage. 
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II .  PROBLEM  FORMULATION 

Consider  Che  geometric  description  of  a  missile  attack  depicted  in 
figure  1.  Initially,  the  weapon  is  launched  from  own  ship  and 

A  A  A 

directed  to  an  aimpoint  specified  by  the  vector  X^  *  (Xw,  Y^). 

However,  because  of  inherent  errors  in  the  missile  delivery  system, 
the  weapon  will  not  necessarily  impact  at  its  intended  aimpoint.  To 
accommodate  possible  discrepancies  between  these  two  points,  the  actual 
weapon  position  at  time  of  detonation  will  be  denoted  by  *  (Xw>  • 
For  computational  convenience,  all  errors  associated  with  weapon 
positioning  are  tacitly  assumed  to  be  Gaussian-distributed  with  the 
following  statistics: 


E{X  } 


(la) 


E{ (X  -  X  )  (X  -  X  )'} 
-v  — w  -v  -v 


w 


(lb) 


where  the  prime  symbol  '  is  used  to  denote  matrix  transposition,  and 
E  {*}  is  the  statistical  expectation  operator. 
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Figure  1.  Geometric  Description  of  a  Missile  Attack 
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In  a  similar  manner,  the  actual  and  estimated  values  of  target 
position  at  time  of  weapon  detonation  will  be  defined  by  X^  =  (Xt>  Y  ) 

A  A  A 

and  X^  *  (Xfc,  Yt),  respectively.  To  determine  the  statistical 
characteristics  of  Xt>  observe  (in  figure  2)  that 


X„  =  X_  +  Xk 
— t  —to  — ta 


(2) 


where  X  represents  target  position  prior  to  alertment  of  a  missile 

attack.  Estimates  of  the  vector  X  are  typically  obtained  from 

shipboard  target  motion  analysis  algorithms  which  may  also  provide 

statistical  descriptions  of  the  associated  estimation  errors.  Since 

these  errors  are  modelled  as  Gaussian  random  variables  (reference  1), 

X^  will  also  be  Gaussian-distributed  with  mean  and  covariance  given 
—to 

by 


E'2to>  -  ho 


(3a) 


*t(i»  '  (Xbo  '  "  ho 


(3b) 
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Figure  2.  Vector  Diagram  of  Relative  Target  Position  at  Time  of 

Weapon  Detonation 
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The  vector  X  defines  the  target  position  at  time  of  weapon 
detonation  relative  to  its  position  at  time  of  alertment. 
Unfortunately,  errors  associated  with  estimating  X  have  a 
complicated  non-Gaussian  probability  distribution  that  functionally 
depends  upon  the  target  evasion  characteristics.  To  devise  tractable 
approximations  for  this  distribution  requires  a  suitable  geometric 
model  of  the  target  escape  region.  The  simplest  model  is  a  circle 
centered  at  whose  radius  increases  with  time  in  proportion  to 

the  maximum  target  speed.  Although  such  models  are  frequently  used  in 
practical  applications,  they  are  somewhat  unrealistic  because  they 
implicitly  assume  that  the  target  has  a  zero  turning  radius. 

Recently,  more  realistic  techniques  of  modeling  and  target  evasion 
characteristics  have  been  developed.  The  resulting  escape  regions, 
shown  in  figure  3,  accurately  reflect  constraints  on  the  target 
dynamics  by  utilizing  a  non-zero  turning  radius  in  conjunction  with  a 
finite  maximum  speed.  At  any  prescribed  time  after  alertment,  the 
target  must  reside  within  a  region  whose  boundary  consists  of  an 
involute  (dotted  curve),  and  arcs  of  the  two  minimum-turning-radius 
circles . 
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As  an  example,  the  darkest  shaded  area  in  figure  3  depicts  the 
target  escape  region  1  minute  after  alertment.  In  this  figure,  the 
origin  of  coordinates  is  defined  by  the  vector  X  (the  assumed 
target  position  before  alertment). 

Because  the  actual  escape  region  is  difficult  to  describe 

analytically,  it  will  be  approximated  by  a  circumscribing  circle  as 

shown  in  figure  4.  Although  this  approximation  introduces  errors  into 

subsequent  probability  computations,  these  errors  decrease  with  time 

inasmuch  as  the  actual  and  approximate  escape  regions  eventually 

coincide  (reference  2).  For  convenience  of  notation,  R  is  used  to 

Q 

denote  the  radius  of  the  circumscribing  circle,  and  X  defines  the 

L  3 

location  of  its  center  relative  to  X_  . 


Figure  4.  Circular  Approximation  of  Target  Escape  Region 
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While  it  is  known  that  an  alerted  target  must  necessarily  reside 
within  a  prescribed  escape  region,  more  specific  data  to  further 
refine  the  localization  process  is  usually  unavailable.  Evasive 
maneuvers  are  almost  certain  to  occur  because  of  the  crisis  nature  of 
the  situation  (i.e.,  the  target  knows  it  is  under  attack);  however, 
the  absence  of  additional  information  precludes  classifying  one  type 
of  maneuver  as  more  likely  thav  another.  In  light  of  this,  it  is  not 
unreasonable  to  assume  that  the  errors  associated  with  estimating 
X  are  uniformly  distributed  over  the  prescribed  escape  region. 

Under  such  an  assumption,  it  can  be  shown  that 


(4a) 


Et(2ta 


*«>•>  ■  <V2,Jl 


(4b) 


where  I  is  the  two-dimensional  identity  matrix,  and  the  pertinent 
statistical  computations  are  performed  using  a  circular  approximation 
as  depicted  in  figure  4. 

To  determine  the  probability  of  target  damage  resulting  from  a 
missile  attack,  it  is  necessary  to  mathematically  combine  the 
statistical  information  embodied  in  equation  sets  (1)  through  (4)  with 
known  weapon  damage  characteristics.  This  task  may  be  accomplished 
via  application  of  the  law  of  total  probability  (reference  1),  which 
leads  to  the  expression 

Prob[target  damage]  ■  Prob[target  damage jR]p(R)dR  (5) 


where 


R  *  |^|  *  range  from  target  to  weapon  at  time  of  detonation 
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ft  *  set  of  all  possible  values  of  R. 

Here,  the  conditional  probability  function,  Prob  [target 
damage | R],  defines  the  probability  of  target  damage  given  that  the 
relative  range  from  weapon  to  target  is  known  at  time  of  detonation. 

It  is  this  term  that  incorporates  weapon  damage  characteristics  into 
the  effectiveness  prediction  model.  The  remaining  term,  p(R), 
represents  the  probability  density  of  R  and  is  used  to  account  for 
target  localization  errors,  target  evasion  character isitcs ,  and  weapon 
aimpoint  delivery  errors.  To  see  this,  observe  from  figure  2  that  the 
vector  X  can  be  written  as 

X  =  X„  +  X„  -  X  .  (6) 

—  —to  — ta  — w 

Thus,  both  the  magnitude  and  direction  of  X  will  be  affected  by  target 
localization  errors  via  XtQ,  target  evasion  characteristics  via 
Xfca,  and  weapon  placement  errors  via  X^. 

To  compute  the  probability  of  target  damage  on-line  and  in  real 
time,  it  is  first  necessary  to  explicitly  represent  p(R)  and 
Prob  [target  damage  Jr]  as  integrable  functions  of  R.  Unfortunately, 
it  is  virtually  impossible  to  develop  exact  mathematical 
representations  for  these  two  functions  which  are  also  amenable  to 
closed-form  integration.  Nevertheless,  suitable  approximations  can  be 
employed,  thus  allowing  equation  (5)  to  be  evaluated  without  recourse 
to  numerical  integration.  The  desired  approximations  for  p(R)  and 
Prob [target  damage | R]  are  computed  in  the  following  section. 
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III.  APPROXIMATION  OF  THE  PROBABILITY  DISTRIBUTIONS 
DETERMINATION  OF  THE  PROBABILITY  DENSITY,  p(R) 

To  mathematically  describe  the  probability  density,  p(R),  it  is 
convenient  to  utilize  the  concept  of  characteristic  functions 
(reference  1).  In  particular,  if  p(x)  defines  the  probability 
density  of  a  two-dimensional  random  vector  x'  *  (x,  y),  then  this 
density  and  its  associated  characteristic  function  $x(u)  satisfy 
the  relations 


p(x) 


1 

(2ir)2 


(7a) 


*  (<u) 

X  — 


f  f  p(x)e&'2dx 
>y  —co  J 


(7b) 


where  OJ  is  a  two-dimensional  vector  of  the  form  cu' 


(ojx,  OJ  ). 


For  the  problem  under  consideration,  observe  from  equation  (6) 
that  the  vector  X  is  a  linear  function  of  X_  ,  X„  ,  and  X  . 

Since  the  errors  associated  with  estimating  these  latter  three  vectors 
are  uncorrelated,  it  follows  (from  reference  1)  that  the 
characteristic  function  of  p(X)  statisfies  the  relation 


4  (w)  -  4  (w)4  (a>)4  (-»)  (8) 

a  —  to  —  ta  —  w  — 

where  4  (w),  4_  (w),  and  4  (u>)  are  characteristic  functions  affiliated 

to  —  ta  —  w  — 

with  the  probability  densities  p(XtQ),  p(Xtft),  and  p(X^) , 
respectively.  In  view  of  the  assumption  that  Xfc0  and  X  are 
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Gaussian-distributed ,  the  functions  $  (w)  and  $  (w)  take 

to  —  w  — 

the  form  (from  reference  1) 

■  exp(jaj'Xt- 
~  1 

$  (uj)  ■  exp(joj'X  -  rto'P  w)  . 
w  —  r  J - w  2—  w— 


(9) 

(10) 


As  discussed  previously,  the  vector  X  is  presumed  to  be 
uniformly  distributed  over  the  target  escape  region;  thus, 


p«t.)  • 

Although  this  expression  may  be  utilized  without  alteration,  it  will 
facilitate  subsequent  computations  to  approximate  p(Xta)  via  the 
formula 


R 


X^  -  X  I  <  R 
ta  — ta*  —  e 


(11) 


X„  -  X  >  R 
i— ta  — ta1  e 


p„<2 


ta 


2 

exp{-(n+l)p  } 

^  (n+l)k 

UJ 

PR2 

e 

^  * 

.2k 


(12a) 


Here,  n  determines  the  order  of  the  approximating  series,  and 


A 


e 


(12b) 


The  development  of  the  finite-sum  Gaussian  distributions  described 
in  equation  set  (12)  is  derived  in  appendix  A.  They  are  often  used 


10 


TS  5822 


to  evaluate  complicated  probability  integrals  in  terms  of  elementary 
functions.  It  might  be  noted  that  application  of  equation  set  (12), 
in  lieu  of  equation  (11),  does  not  change  the  mean  and  covariance  of 
Xta»  Consequently,  equation  set  (4)  will  continue  to  provide  a 
valid  statistical  description  of  X  regardless  of  which  density 

*"“Ca 

function  is  chosen  to  approximate  the  pertinent  error  distribution. 
For  the  case  n=l,  observe  that  p^(Xta)  is  Gaussian  and  reduces  to 


(13) 


Finally,  the  limiting  approximation  as  n  approaches  infinity  is  the 
original  uniform  distribution;  i.e., 


lim  ^n(— ta)  "  p(4a)  *  <14) 

n-*» 


The  characteristic  function  for  p  (X„  ),  developed  in  appendix 

n  — ta 

B,  is  given  by 

<&ta^— ;n^  =  n  [lpeu,2]exp^— ^ '-ta  ~  \  (15a) 


where 


U) 


j  to 


V 


2^2 
U)  +  uj 
x  y 


(15b) 


R2 

e 

2(n+l) 


(15c) 


and  are  Laguerre  polynomials  (from  reference  2),  defined  by 
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^P)C0-  t  (pij  h  <-5>k  a. 

Here,  n  appears  as  an  argument  in  the  characteristic  function  of 
Pn(Xta)  to  denote  the  order  of  the  approximating  series. 

Since  equations  (8),  (9),  (10),  and  (15a)  completely  specify  the 
characteristic  function  of  p(X),  this  probability  density  may  be 
determined  analytically  from  equation  (7a)  as  follows: 


p(X) 


1 

(2n)2 


J~  J'  $x(a>)e  ■*—  ^dw  . 

CO  00 


(17 


To  obtain  the  distribution  of  range  errors,  recall  that  R=)x] . 
Consequently,  if  the  cartesian  coordinates  of  X  are  given  by  X'  = 
(X, Y),  the  associated  polar  coordinates  may  be  expressed  as 


R  = 


Vx2  +  Y2 


(18a 


0  =  tan"1  (Y/X) 


where  0  <  0  <  2u 


(18b 


The  Jacobian  of  the  transformation  from  cartesian  to  polar  coordinates 
then  takes  the  form 


3X 

3Y 

3R 

3R 

COS0 

-Rsin0 

J  =■  det 

3X 

3Y 

=  det 

sin0 

Rcos0 

30 

30 

- 

- 

Using  the  transformation  theorem  for  probability  density  functions,  it 
can  be  shown  (from  reference  3)  that 
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DETERMINATION  OF  CONDITIONAL  TARGET  DAMAGE  PROBABILITY, 

Prob [target  damage |r] 

To  mathematically  describe  the  probability  function  Prob  [target 
damage) R]  requires  explicit  knowledge  of  the  weapon  damage 
characteristics.  For  this  study,  missiles  with  circular  interdiction 
zones  are  used  as  an  example  to  determine  the  probability  of  target 
damage.  These  characteristics  will  be  specified  by  concentric 
circular  regions  as  shown  in  figure  5  and  defined  as  follows: 

1.  Lethality  or  Kill  (Seaworthiness  Impairment).  Lethality  from 
underwater  explosions  is  associated  with  overpressure  necessary  to 
rupture  the  submarine  pressure  hull. 

2.  Immobilization  (Mobility  Impairment).  Immobilization  is 
associated  with  shock  damage  to  internal  equipment.  It  is  defined  as 
that  condition  whereby  the  target  is  rendered  virtually  immobile 
because  of  damage  to  propulsion  machinery. 

3.  Mission  Abort  (Weapon  Delivery  Impairment).  Weapon  delivery 
impairment  is  associated  with  shock  damage  to  internal  equipment 
involved  in  target  detection  and  weapon  delivery.  It  is  defined  as 
that  condition  whereby  the  target  is  virtually  unable  to  detect  enemy 
vessels  or  release  weapons  effectively. 

Each  region  is  associated  with  a  50  percent  probability  of  target 

damage,  and  the  respective  centers  all  coincide  with  weapon  position 

at  time  of  detonation.  Note  also  that  regions  of  lesser  potential 

destruction  always  encompass  regions  of  greater  potential  destruction 

as  subsets.  For  computational  convenience,  R,  will  be  used  to 

a 

parametrically  define  the  radius  of  the  pertinent  damage  envelope 
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under  consideration.  In  this  way,  the  probability  of  the  missile 
inflicting  any  prescribed  type  of  damage  to  a  target  can  be  determined 
from  a  single  formula  by  simply  assigning  the  appropriate  value  to 


Figure  5.  Geometric  Description  of  Missile  Damage  Characteristics 

Although  specific  numerical  data  are  available,  an  exact 
mathematical  description  of  Prob[target  damage|R]  does  not  exist. 
However,  from  practical  considerations  it  is  apparent  that  this 
function  will  decrease  monotonically  as  R  increases,  and  also  must 
satisfy  the  inequality 

0  <  Prob[target  damage) R]  <  1  (23) 

for  0  <  R  <  ®.  These  characteristics  suggest  that  finite-sum 
Gaussian  distributions  again  be  used  as  a  basis  for  approximating 
Prob [target  damage |r];  i.e., 


Prob[target  damage|R]  *■  Prob[target  damage|R,m] 


(24) 


where 


Prob [target  damage j R,m]  =•  exp(-m^2)  ]T  ry-  52k.  (25) 

k-0 

As  before,  m  determines  the  order  of  the  approximating  series,  and 

5  =*  R/Rd  •  (26) 


Graphs  of  equation  (25)  for  m-1,2,3, . . .10  and  m=®  are  plotted 
in  figure  6.  As  expected,  these  functions  all  decrease  monotonically 
as  R  increases.  Note  also  that  the  limiting  conditional  probability 
may  be  expressed  in  the  form 


Prob[target  damage |R,«] 


1  0  <  R  <  R , 

—  —  a 

0  8  >  Bd 


(27) 


Because  of  its  simplicity,  equation  (27)  is  frequently  used  in 
practical  applications  even  though  it  is  somewhat  artificial.  Indeed, 
such  an  approximation  implicitly  assumes  that  missile  attacks  will  have 
uniformly  devastating  effects  over  the  entire  weapon  damage  envelope. 

In  light  of  this,  it  is  perhaps  more  realistic,  and  certainly  more 
economical,  to  employ  lower-order  approximations  when  possible. 
Accordingly,  the  parameter  m  will  not  be  specified  here  so  as  to 
enhance  modelling  flexibility  and  realism. 


I 
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NORMALIZED  RANGE  FROM  WEAPON  TO  TARGET  (R/Rd) 


Figure  6.  Plots  of  the  Function  -  Prob£ target  damage]R,mJ 
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IV.  COMPUTATION  OF  THE  TARGET  DAMAGE  PROBABILITY 

The  probability  of  target  damage  resulting  from  a  missile  attack 
may  now  be  determined  analytically  by  utilizing  approximations 
developed  in  section  III  in  conjunction  with  the  law  of  total 
probability.  More  precisely,  substituting  equations  (22),  (24),  (25), 
and  (26)  into  equation  (5),  and  evaluating  the  resultant  integral  over 
ft,  (see  appendix  C)  leads  to  the  expression 


00  00 

Prob  [target  damage]  *  If  $x(u)'f(aj)d w 


where 


^  -  4 


[l »/] 


»d  ■  Ed/2" 


(29a) 


(29b) 


and  L^_^[£]  is  the  Laguerre  polynomial  defined  previously.  In  deriving 
equation  (28),  it  might  be  noted  that  integration  of  the  exact 
probability  distributions  over  ft  is  mathematically  equivalent  to 
integrating  the  approximate  distributions  over  a  semi-infinite 
interval,  0  <  R  <  «. 


If  equations  (8),  (9),  (10),  (15a),  and  (29a)  are  combined,  the 
integrand  in  equation  (28)  may  be  rewritten  as 

$x(wMu>)  «  ^  Ln-l[lPew2]Lm-l[Fd^2]exp 
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t 

► 
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t 

i 

i 


where 

Ola) 


(31b) 

and  I  is  a  (2x2)  identity  matrix.  Utilizing  equation  (16),  in 
conjunction  with  equations  (28)  and  (30),  then  yields 

Prob [target  damage]  *  ^  pj^+l  (32) 

where 

\*i  =  4r  'T7TT  Jj®  -  kpx^  •  (33) 


X  = 


x,  +  x„  -  x 

—to  — ta  — w 


P  = 
x 


Px(l,l)  Px(l,2) 
Px(l,2)  Px(2,2) 


*  P„„  +  P  +  (p  +  p.)I 
to  w  e  d 


The  last  remaining  task  is  to  determine  ftj^  in  terms  of  the 
pertinent  modeling  parameters.  This  may  be  accomplished  by  employing 
equation  (15b)  and  the  Binomial  Expansion  Theorem  to  explicitly 
evaluate  the  integral  appearing  in  equation  (33)..  Pertinent 
mathematical  details  are  presented  in  appendix  D,  and  the  final  result 
takes  the  form 


ft 


k+i 


exp{-5s(C2  +  0} 

_ 2 _ i _ 

Vx+  x_ 


k+i 


(34) 
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1 .  Target  Localization  Data 


estimated  target  position  prior  to  alertment 
of  a  missile  attack;  this  vector  is  defined 
with  own  ship  as  the  origin  of  coordinates 

PT0(1 ,1)  PT0(1,2) 

PTO (2,1)  PTO(2,2) 

estimated  target  position  after  alertment  to 
a  missile  attack;  this  vector  is  defined  with 

X^  as  the  origin  of  coordinates 
to 


RE  *  Rfi  evasion  radius  of  target 


integer  which  determines  the  order  of  the  Gaussian 
sum  approximation  to  p(Xta)  / 


'XTA(l)' 

XTA(2) 
u  - 


=  P 


to 


covariance  matrix  of  X_ 

—to 


XT0(1) 

XT0(2) 


2.  Weapon  Placement  Data 


XW(1) 

_XW(2) 


estimated  weapon  position  at  time  of  detonation; 
this  vector  is  defined  with  own  ship  as  the 
origin  of  coordinates 


PW(1,1)  PW(1,2) 

PW(2, 1)  PW(2,2) 


covariance  matrix  of  X 

— w 
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3.  Weapon  Damage  Characteristics 

RD  *  R^  damage  radius  of  weapon 


s  integer  which  determines  the  order  of  the  Gaussian 
m  sum  approximation  to  Prob[ target  damage  Jr] 

Once  the  preceding  data  have  been  specified,  the  algorithm  will 
compute  PROB  »  probability  of  target  damage  each  time  the  subroutine 
missile  (PROB)  is  called. 
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V.  SUMMARY 

A  real-time  procedure  for  predicting  the  effectiveness  of  a 
missile  attack  against  prescribed  target  threats  has  been  developed. 
The  required  probability  of  target  damage  is  specified  in  terms  of 
known  or  measurable  weapon/ target  characteristics  and  associated  error 
statistics.  Furthermore,  the  computational  algorithm  is  sufficiently 
compact  so  as  to  allow  weapon  performance  predictions  to  be  updated 
on-line  as  tactical  conditions  change.  As  a  result,  the  algorithm  has 
potential  value  as  an  automatic  decision  aid  to  assist  in  the  optimal 
deployment  of  missiles. 

Although  the  weapon  damage  characteristics  and  target 
localization  solution  are  constant  inputs  to  the  probability  model, 
other  inputs  (e.g.,  missile  aimpoint  coordinates  and  target  evasion 
characteristics)  are  time-dependent  functions  which  change  as  the 
missile  aimpoint  and  flight  time  are  varied.  In  view  of  these 
properties,  it  is  evident  that  an  optimum  aimpoint  exists  which  will 
maximize  targt  damage  probability  for  any  prescribed  tactical 
scenario.  Future  work  will  focus  on  application  of  stochastic 
optimization  techniques  to  solve  the  non-linear  equations  which 
determine  the  optimum  missile  aimpoint. 
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APPENDIX  A 


DEVELOPMENT  OF  FINITE-SUM  GAUSSIAN  DISTRIBUTIONS 


Given  a  uniform  density  of  the  form 


7TR 


X_  \  <  R 

ta  — ta1  —  e 


I— ta  ~Xta‘  >Re 


(A-la) 


where  X_ 
— ta 


(A-lb) 


and  X„  »  E 
— ta 


N 


(A-lc) 


it  is  desired  to  approximate  it  by  a  finite-sum  of  Gaussian 
distributions  of  the  form 


P«t„>  =  •  *  «P 


ta 


-ta1 


n-1 


1*0 


a1 lx.  -  X„ 
'-ta  — ta 

i! 


2i 


(A-2) 


In  the  approximating  series,  n  determines  the  order  of  the  series.  As 
n  approaches  infinity 


P„<5ta)  “  p(2ta)- 


ta 


(A-3) 


The  constants,  K  and  a  in  equation  (A-2)  are  found  by  evaluating  the 


A-l 
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statistics  of  the  approximating  density  so  that  the  mean  and 
covariance  between  the  two  densities  remain  unchanged.  Also,  for 

?n ta^  t0  a  probability  density  it  must  satisfy  the  following 

constraints : 


)  >  0 


(A-4a) 


and 


<*> 

J 


)  dX„ 


1. 


(A-4b) 


Before  constraint  (A-4a)  can  be  satisfied,  the  constants  K  and  a 
must  be  found.  First,  it  will  be  shown  that  the  approximating 
probability  density  satisfies  constraint  (A-4b).  In  order  to  simplify 
the  integration  involved  in  satisfying  the  constraint,  a 
transformation  from  rectangular  to  polar  coc  dinates  is  necessary. 
Using  figure  A-l  let 

A 

x-x  ■  R  cos  0 


and 

A 

y-y  =  R  sin  0 


or 

R2  -  (x-x)2  ♦  (y-J)2 


and 

0  ■  tan 


(y-y)/(x-x) 


A- 2 


Figure  A-l.  Relationship  Between  Polar  and  Rectangular  Coordinates  of 
the  Actual  and  Estimated  Target  Position 


Substituting  for  x  and  y  into  equation  (A-2)  yields 


vr  *  \  >/',  x  .  _  -otR 

Pn(Xta)  -  PQ(x,y)  *=  Pn(R,0)  *  Ke 


2  aS2i 


(A-5) 


1*0 


From  reference  A-l,  the  integration  in  polar  coordinates  becomes 


/  /  ■ K  £  fr  /  f 


“-1  -1  ^  -oR22i 

6  K 


9(x,y) 
3  dR,0) 


dRd©  *  1 


where  the  Jacobian  of  the  transformation  is  defined  by 
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J 


3(x,y) 


a(R,e) 


det 


“3* 

3R 

3x 

.30 


h. 

3R 


=  R. 


Integrating  over  0,  the  equation  reduces  to 


2ttK 


f 


-OR 


R2l+1  dR 


1. 


1=0 


(A-6) 


From  reference  A-2  (p.65,  #314,  2b)  the  integral  in  equation  (A-6) 
becomes 


2irK 


E 


n-l 


i=o 


(A-7) 


The  summation  in  equation  (A-7)  can  be  written  as 


n-l 

y*  1  =  l+l+l+l . 1  *  n.  (A-8) 

i”o 


Substituting  equation  (A-8)  into  equation  (A-7)  and  solving  for  K 
yields 


K  *  Qi/Tin. 


(A-9) 


Selecting  the  constant  K  by  equation  (A-9)  insures  that  the  constraint 
given  by  (A-4b)  is  satisfied. 
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Next  the  mean  of  the  approximating  density  is  found.  The  mean  is 
given  by 


mean 


00  00 

/  /  X_  p  (X_  )  dX 

J  J  — ta  *n  — ta  -ta 


which  can  be  written  as 


mean 


09  00 

/  feu'  i 


00  00 

ffi 


p  (X,.  )dX  .  (A- 
ta  rn  — ta  — ta 


-00  -CO 


Utilizing  constraint  (A-4a) ,  the  second  term  reduces  to 


X„ 

— ta 


CO  CO 

/A 


(X  )dX  *  X  . 
ta  — ta  — ta 


(A- 11) 


To  evaluate  the  first  integral,  the  substitutions 


x' 

ta 

ta  ta 

dx’ 

»  dx^ 

ta 

ta 

A 

y' 

7ta 

"  yta"  yta 

and 


dA  *  d?ta 


A-! 
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are  made. 


The  integral  now  becomes 


00  03 


If 


CO  00 


If 

—CO  -CO 


ta 


ta 


(x‘  + 

n  ta 


ta*  yta 


yM>d*« 


Making  a  transformation  to  polar  coordinates  where 


x'  ■  R  cos  0 
ta 

y '  =  R  sin  0 

■'ta 


results  in  two  integrals  of  the  form 


2TT  00 

f  f  R  cos  0  Pn(xta  +  xta’  yta  +  RdRdQ 

o  o 


and 


2n  « 

/  /  R  sin  0  p  (x'  +  x 

y  y  n  ta 

o  o 


yta  +  yta)  RdRd0* 
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By  carrying  out  the  actual  integration,  it  is  found  that  both 
integrals  are  zero.  Therefore,  the  mean  of  the  approximating 
probability  density  function  is  the  same  as  the  uniform  probability 
density  as  seen  by  equation  (A-ll). 

The  covariance  of  the  uniform  density  is  given  by 

P„  -  (R2/2)I  (A-15) 

ta  e 

where  I  is  a  2x2  identify  matrix.  The  constants  K  and  ct  in  the 
approximating  density  function  will  be  determined  so  that  its 
covariance  is  also  P  .  The  covariance  is  given  by 

Id 

OO  CO 

covariance-  f  /  (X,.  -  )  (X.  -  X.  )  'pCX.  )dX.  .  (A-16) 

/  /  — ta  — ta  —ta  — ta  n  — ta  — ta 


Making  the  transformation  to  polar  coordinates,  the  result  is 


2ir  ® 


covariance 


J *  J“  [Rcos0  Rsin0]  p^CRcosG,  Rsin0)RdRd0  (A-17) 


o  o 


which  can  be  written  as  four  integrals  of  the  form 


217  /< 

// 


R2cos20  p  (Rcos0,  Rsin0)RdRd0 
u 


(A-18a) 


o  o 
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o 


CO 

2  s 

R  cos©  sin0  p^(RcosG,  Rsin0)RdRd0 


00 

2  wT 

R  sin©  cos©  p^(Rcos©,  Rsin©)RdRd© 


o 


CO 

R^sin^Q  p^CRcos©,  RsinQ)RdRd©  . 


Using  the  trigonometric  identity, 
sinQ  cosG  =  -|sin20, 

the  integrals  given  by  equations  (A-18b)  and  (A-18c)  are  zero, 
integrals  given  by  equation  (A-18a)  and  (A-18d)  can  be  easily 

evaluated  by  substituting  equations  (A-5)  and  (A-9)  for  p(Rcos 
After  the  substitution,  equations  (A-18a)  and  (A-13d)  become 


a 

trn 


n-1 

£ 

i=o 


00  2tt 

7 \Jf' 


e-aR  r2i+3  cos20  dQdR 


A-8 


(A- 18b) 

(A-ISc ) 

(A-18d) 

The 

,  RsinG). 

(A-19a) 


Integrating  with  respect  to  8,  both  equations  (A-19a)  and  <A-19b) 


So- 
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Since  the  covariance  of  the  uniform  density  and  the  approximating 
density  must  be  the  same,  equation  (A-20c)  is  set  equal  to  equation 
(A-15)  (covariance  matrix  of  the  uniform  density),  and  the  constant 
i  is  determined.  The  result  is 

a  =  (n+l)/R2  .  (A-21 ) 

e 

Substituting  equation  (A-21)  into  equation  (A-9),  the  constant  K  is 
given  by 


e 


and  the  approximating  density,  equation  (A-2),  becomes 


exp 


(X„  > 
n  — ta 


n+1 

n 


"(n+1)  lit.'  it. I 


1TR 


n-1 

£ 

i=o 


(n+i)l|x  - X 

ta  — ta 


2i 


i!  R 


2i 


(A-23 ) 


Since  n,  R  ,  and  lx  -X„  I  are  all  greater  than  or  equal  to  zero, 
’  e  '  — ta  — ta 1 

p(Xta)  >0 

and  the  constraint  given  by  equation  (A-4a)  is  satisfied. 
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APPENDIX  B 

DEVELOPMENT  OF  THE  CHARACTERISTIC  FUNCTION 


Given  a  probability  density  of  the  form 


-(n+1)  lx  -  X„  I 
ta  — ta1 


exp 


•r  .  .  n+1 

p  (X^  )  *  - 

n  — ta  n 


R 


irR 


n-1 

E 

i=o 


(n+1)1|^ta^ta 


2i 


i!  R 


2i 


(B-l) 


the  characteristic  function  can  be  found  from  the  equation 


jM'I 


ta  dXta 


(B-2) 


where 


u> 


and  to  * 


Substituting  equation  (B-l)  into  equation  (B-2)  and  letting 


X  »  x„  -  x_ 
—  — ta  — ta 


dX  -  dX„ 
—  — ta 


results  in  a  characteristic  function  of  the  form 
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CO  00 


4>  (to)  =  e-'—  — ta  f  fp  (X  +  X  )  e-^—  -  dX. 
ta  —  J  J r n  -  — ta  — 


-CO  —CO 


Making  a  transformation  to  polar  coordinates,  where 


(3-3) 


0)  COS<j> 

R  cos© 

(0  = 

X  = 

to  sin4> 

R  sin© 

- 

and 

(0  * 

IHI 

R  = 

llxll 

the  characteristic  equation  is  given  by 


$ 


A 


exp  { jSujcos (©~4>)}  RdRd©, 


(3-4) 


where  Rg  is  the  radius  of  the  evasion  circle.  Integrating  with 
respect  to  0  leaves 


(to) 

ta  — 


(n+l)k 

k! 


a 

/ 


exp  ! - ( n+ 1 ) 


(R/R^) 


(R/R  )2k+1  J 
e  o 


^wv)  r-  <5-5> 

e 


B-2 


where  J  iu)R  (R/R  )l  is  the  zeroth  order  Bessel  function  (reference  7). 
o«  e  e  > 


Let  pg  *  R/Rg,  then  equation  (B-5)  becomes 


$  (u>)  =  e^—  -ta  i  £ 

ta  —  n  . 


11-1  (n+l)k 


k=o 


k! 


CO 

/ 


2k 

pe  exp 


-(n+l)pk 


J  (uR  p  )  (n+l)2p  dp  . 
o  e  e  e  e 


(B-6) 


Making  another  substitution,  X  =  (n+l)pg,  equation  (B-6)  becomes 


•«.<«>  -  i  g  4r  / ',(2  ^£) 


dX 


which  integrates  to  (from  reference  A-2,  p.  721,  #6.643,4) 


ta 


„  2D2 

n-1 

’  2p2 

1 

~  a)  R 

iw'X±  e 

YV°> 

oj  R 

a 

-  -  exp 

~  ~ta  4(n+l) 

k=o 

4(n+l) 

(B-7) 


where  are  Laguerre  polynomials  (reference  7).  By  noting  (from 

reference  A-2,  p.  1038,  #8.974,3) 
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n-1 


(o) 


r  2_2  i 

0)  R 

e 

=  L(1) 

r  2  2  i 

uTR 

e 

4 ( n+ 1 ) J 

Ln-1 

_4(n+l) 

equation  (B-7)  becomes 


$  (aj)  =  —  exp 
a  —  n  r 


2_  2  i 

t 

2_  2 

OJ  R 

jo)  x  _  e 

L(1) 

0)  R 

e 

ta  4(n+l) 

1  11 

4(n+l) 

B-4 
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APPENDIX  C 

PROBABILITY  OF  TARGET  DAMAGE 

To  determine  the  probability  of  target  damage,  we  use  the 
expression 

Prob[target  damage]  a  I  Prob[target  damage |r]  p(R)  dr  (C-l) 

or 

where 

R  *  | X |  *  range  from  target  to  weapon  at  time  of  detonation 

*  set  of  all  possible  values  of  R, 

The  conditional  probability,  Prob[target  damage |r],  is  approximated 
by  a  finite-sum  Gaussian  distribution  given  by 

Prob[target  damage|R]  *  Prob[target  damagejR,m] 


(C-2) 


where  m  determines  the  order  of  the  approximating  series  and  R^  is 
the  radius  of  the  pertinent  damage  envelope.  The  probability  density 
function  p(R)  is  written  as 
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p(R) 


_1^_ 

2tt 


00  ao 


If 

-40  -00 


$  (w)R  J  (ojR)  dai 

x  —  o  — 


(C—  3 ) 


where  <J>  (w)  is  the  characteristic  function  defined  by 
x  — 


<t>  (w)  =  4.  (0))$ 

x  —  t  — 

o 


(u)4  (-U)) 
ta  —  w  — 


(C-4) 


and  J  (u)R)  is  the  zeroth  order  Bessel  function  (reference  7). 

0 

Substituting  equations  (C-2)  and  (C-3)  into  equation  (C-l)  yields 


Prob  [target  damage] 


2ir 


m- 1  00  <30 

Ek/ 1''--  P 

k=0  -oo  -oo  Q 


*  J  (oiR)RdRdo) 

o  — 


(C-5) 


Let 


Coj  ) 


m-1  . 

Zb 

k=o 


e-m(R/Rd)  j^(ajR)RdR  ^ 


(C-6a) 


and  carry  out  the  integration.  Since  the  conditional  probability 
distribution  is  an  approximation,  it  is  equivalent  to  integration 
over  the  interval  0  <  R  <  ®  .  Equation  (C-6a)  then  becomes 


C-2 
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^(cu) 


1_ 

2tt 


-m(R/Rd> 


2 

J  (uiR)RdR  . 
o 


(C-6b) 


The  integral  in  equation  (C-6b)  has  the  same  form  as  equation  (B-6). 
By  using  the  substitution, 


X  *  m(R/Rd)2 


the  solution  to  equation  (C-6b)  is  given  by 


iKw) 


m-L 


R 


d 

2m 


exp 


2  2 
-a)  Rd 


4m 


(C-7) 


Using  equation  (B-8),  equation  (C-7)  can  be  written  as 


1_ 

2ir 


2„  2  ' 

2 

2_2 

“  Rd 

Rd 

0)  Rd 

4m 

2m  exp 

4m 

.. 

(C-8) 


Substituting  equation  (C-8)  into  equation  (C-5),  the  probability  of 
target  damage  now  takes  the  form 


00  00 

Prob  [target  damage]  ■  J "  ^$x(u)  ty(u)du. 


(C-9) 
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APPENDIX  D 


EVALUATION  OF  Q 


Determination  of  the  probability  of  target  damage  requires 
explicit  evaluation  of  the  following  integral: 


Vi-ir  tHt-  /  f  [rp“p  <I>-1) 


where 


X_  +  X_  -  X  , 
—to  — ta  -w 


P^  is  a  symmetric  positive  definite  covariance  matrix  of  the  form 


Px(I,l)  Px(l,2) 


Px(l,2)  Px(2,2) 


Pto  +  Pw  +  (pe  +  pd)l» 


and  I  is  a  (2x2)  identity  matrix.  A  diagonal  matrix  which  contains 
the  eigenvalues  of  Px  can  be  formed  by  using  the  formula 


A  -  M  *-PxM  -  diag  (Xx,Ay) 


(D-2) 
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where  M  is  the  modal  matrix  and  X  , 

x’ 

Solving  equation  CD— 2 )  for  P^  yields 


X 

y 


are  the  eigenvalues  of  P  . 


P 


x 


MAM 


(D-3) 


Since  P^  is  symmetric,  the  eigenvectors  are  orthogonal  ar-’  H  ^  =  M' . 

Now  equation  (D-3)  becomes 

Px  =  MAM' .  (D-4) 

Substituting  N'  =  M  and  Q'Q  *  A  into  equation  (D-4)  results  in 

Px  *  N’Q’QN.  (D-5) 

Substituting  equation  (D-5)  into  equation  (D-l),  the  exponent  becomes 

~  i 

exp  { jtu'X  -  ^  (QNu)'  (QNw)}  .  (D-6) 


Making  another  substitution, 

£  =  QNW  (D-7)  ; 

and  using  equations  (D-5),  (D-6)  and  (D-7),  the  integral  in  equation  ;  ' 

(D-l)  reduces  to  S  . 

c 

II 

D-2 


The  partial  derivatives  take  the  form 
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3u) 

3u) 

_x  .  _ 

X 

3£  xx 

X 

y 

3u> 

3u) 

— 2.  =  T 

35x  yx 

3C 

y 

and  the  Jacobian  of  the  transformation  is  given  by 

3(t o  ,  0)  ) 
x  y 

3<V  O 


T  T 

xx  xy 


T  T 

yx  yy 


T  = 


¥[  [q  " 


l.j 


(D-9) 


The  determinants  in  equation  (D-9)  are  given  by 
iQl  *  VTT  and  |N)  -  l. 


The  second  determinant  was  defined  to  be  the  transpose  of  the  modal 
matrix,  M.  Since  M  ^  *  M’ ,  the  determinant  of  M  must  be  unity; 
therefore,  equation  (D-9)  becomes 


3(u>  ,u>  ) 

x’  y 


VTT 

x  y 


(D-10) 


Substituting  equation  (D-10)  into  equation  (D-8)  and  letting  M  *  N' 
yields 


R 
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2 

y 
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Equation  (D-ll)  is  then  written  as 
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exp 
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(D-12) 


Using  the  Binomial  Expansion  Theorem  in  equation  (0-12)  yields 


11) 


D-5 
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Equation  (D-13)  can  be  written  as 


(-1) 
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2ir  X  k!i! 
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(D-14a) 


where 


k+i-£ 


p  12 

e^x^x  "  2^x  d£ 


(D-14b) 


and 


f,(5  ,  X  ) 

A  *y  y 


~  12 
ejCyCy  -  2  5y  ^ 


(D-14c) 


Making  the  substitution,  £x  *  /2Z  in  equation  (D-14b)  yields 
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Equation  (D-15)  can  be  integrated  to  yield  (from  reference  7) 
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X  1 

where  l  - J  is  a  Hermite  polynomial  given  by 
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Substituting  equation  (D-16b)  into  equation  (D-16a)  and  reversing  the 
summation  yields 


f  (£_,  X  )  «  Viir  exp 
m  x  x 
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The  factorials  (2k)!  and  (2m)!  can  be  written  as 
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then 


fm(5x>  ^x)  =■  /2i"  m!(-l)m  exp 


0  (5  ,  X  ) 

m  ’x  x 


(D-20) 


Solving  equation  (D-14c)  for  f.(£  ,A  )  results  in  a  similar 

jc  y  y 


equation  of  the  form 


wv 


(D-21) 


Substituting  equations  (D-20)  and  (D-21)  into  equation  (D-14a)  yields 


;+i 


yTT 

x  y 


(?) 


k+i 


'k+i\  V7  ~-(£2  +  E2) 

kj  2j  e  2  x  vWWi(VV’  (D- 

4=0 


22) 


D-9/D-10 
Reverse  Blank 


TR 


APPENDIX  E 


FORTRAN  PROGRAM  (WITH  FLOWCHART) 

TO  EVALUATE  THE  PROBABILITY  OF  TARGET  DAMAGE 
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INITIALIZATION 

TARGET  DATA;  £»o,pto12to,RE,  N 
WEAPON  PLACEMENT  OATA;  £W,PW 
WEAPON  DAMAGE  CHARACTERISTICS;  RD,  M 

zzzzm 

COMPUTE  RANGE  BETWEEN  ESTIMATED 
TARGET  AND  WEAPON  POSITION  AT 

DETONATION—  £.  (EQUATION  (31  a)) 
. . . . . " 

COMPUTE  COVARIANCE  MATRIX  ASSOCIATED 
WITH  X  (EQUATION  (3ib>) 

_ : . . l _ _ _ 

COMPUTE  EIGENVALUES  OF  RESULTANT 
COVARIANCE  MATRIX  (EQUATION  (36)) 
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COMPUTE  EQUATION  (39) 
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